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ON THE AUTO AND CROSSCORRELATION 
OF PN SEQUENCES 

i 

James C. Morakis 
ABSTRACT 

The autocorrelation and crosscorrelation proper- 
ties of pseudorandom (PN) sequences are analyzed by 
using some important properties of PN sequences. 
These properties make this discuss ion understandable 
without the need of linear algebraic approach. The 

analysis is followed by some experimental -results . 

! 
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ON THE AUTO AND CROSSCORRELATION 
OF PN SEQUENCES — 


Due to increasing interest in PN codes for applications in multiple user and 
or combat of multipath and others the author has decided to compile some of the 
already known and other recent developments on the correlation properties of 
PN coding schemes and their intercomparison. The basic linear PN code gen- 
eration is discussed extensively in[i] and some of the important properties of 
PN codes will be repeated for convenience. , 


A maximal linear PN Sequence is defined as a sequence generated by a 
linear* feedback shift register m stages long such that the length of the sequence 
is the maximum possible for a constant m, It can be easily shown that if such 
a maximal sequence exists itp'length must be 


/ L = 2" - I 


Figure 1 demonstrates such a feedback shift register with h s either 0 or 1 
h. = 1 means that there is a connection and h. = 0 means that there is no con- 
nection (h tf ~ h m =1 always). One can consider the contents of the shift register 
as the state of the shift register. If the state is known at time t . then the next 
state is given by a relationship of the original state and the transition matrix 
A**. 

The characteristic equation of the matrix A turns out to be E 1 1 


v 


P A (*) 


m 

- Z>. 


» A - y wi th h. = 0 , 1 and h 0 = = 1 • 

vssas i = 0 


■ "f 


* The only allowable logic is mod 2 which Is on exculusive-or gate. 

** See Apendix A for the exact form of A In term* of the connection*. * 
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Figure 1— A Linear Feedback Shift Register 


Tliis polynomial is also called the generating polynomial of the sequence. If 
this polynomial is primitive then the sequence is maximal. This polynomial, 
just like any other polynomial of degree m, has m roots which might be desig- 
nated as a , a j , a a m - 1 * The knowledge of any one of the above m 

roots is sufficient because all the others can be represented as powers of the 
one chosen. Thus we can say that P A (x) has a root a. If a polynomial has a 
root /? such that fi = a k , k any integer, then that polynomial is either polynomial 
' P A (x) or another primitive polynomial. The number of distinct maximal se- 
quences that can be generated by an m-stage shift register is given by [1] 


<£( 2 * - 1 ) 
m 


where 4> (j) is the Euler phi function and it is equal to the number of integers 
prime to j but less than j. So if j is a prime <£ (j) = j - 1. Table 1 gives the num- 
ber of maximal sequences the length L and other information for each m. Table 
1 also gives a listing of the first few Mersenne primes which are defined as the 
m’s which are prime such that 2 m - 1 is also a prime. Obviously Mersenne 
primes result in large numbers of maximal sequences. To investigate the cor- 
relation properties of shift register generated sequences we must first discuss 
some pertinent properties. A sequence of period h can be viewed as a sequence 
of L states (or L n-dimensional vectors; h = 2" - 1 for binary). 

Property A. A maximal binary sequence will occupy all 2 n ~ 1 non-zero states 
before it starts repeating. 


Property B. (Balance property) The number of l’s in a maximal sequence ex- 
ceeds the number of 0's by 1. 


Thus the number of zeros of a maximal sequence is equal to 2 1 - 1 and the 

number of i's is equal to 2 n “ 1 
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Table 1 



m 

L = 2° - 1 

Prime - 
factors 
of 2“ - 1 

2 

3 

3 

• 

3 

7 

■ V 7 ’ 

4 

15 

3*5 

/ . ■ 

5 

31 

31 

. 6 

63 

i '• 3 2 • 7 

7 

127 

\ ’ 127 

8 

- 1 

255 

. \ 3*5*17 

\ 

9 

? 511 

7 • 73 

10 

1023 

3 • 11 • 31 

11 

2047 

23 * 89 

12 

4095 

3 2 • 5 • 7 • 13 

13 

8191 

8191 

14 

16383 ~ 

3 • 43 * 127 

15 

\ 32767 1 

7*31*151 

16 

65535 

5 « 3 • 17 • 257 

17 

131071 

131071 

18 

262143 

3 3 • 7 • 19 • 73 

19 

524287 

524287 

20 

1048575 

3 • 5 2 • 11 • 31 * 41 



(2 - 1 ) 


2 

6 

8 

30 

36 

126 

128 

432 

600 v 


1936 

1728 

8190 

10584 

27000 

32768 

131070 

139968 

524286 

480000 


6 

18 

16 

48 

60 

' 

176| 

144 1 

I 

630 

756 

j 

1800 

2048 

7710 

7776 

27594 

24000 


Mersenne Primes: 

2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 3217, 
4253, 4423, 9689, 9941, 11213, 








Property C. (The Shift -and -add property.) 


... , 


When a maximal sequence is first shifted by t bits and then added to itself ; 
the resulting sequence is the same sequence delayed by some time t* i»e. '■ ; 


?»-■ 


S A (t) + s A ( t - k) = s A ( t Ht,) for k ^ 0 


Property D. This is really a transformation process that applies to all binary 
waveforms. It facilitates the calculation of correlation by transforming multi- 
plication into mod 2 addition. As a result of this transformation l’s in a binary 
sequence are transformed into 0's, the -l’s are transformed into l’s and the 



. V t 


( 1 ) 


where S c is the sequence formed by adding S A to S B mod 2 f bit by bit. 

\ • ■ * ' ' 

The above four concepts will enable us to arrive at some of the most im- 
portant conclusions on the values of auto and crosscorrelation of linear PN 
sequences without becoming dazed by the complex maze of modern math which 
would normally be necessary for this purpose. • 

* **- 1 '* «■ y \ . 

Auto-Correlation of a Linear Maximal PN Sequence ; 

• • ■ j 

The goal is to evaluate ' ;:V;;Ks i; K i 


: ;* • j ? « ; ; . - * , a 's ; ■ * J 


\ Of)*i 2> - k >. ,, 

• " , j - I ’ • : ‘ ' ‘ “ 


where the dimension of t and k is time in integer numbers of bits. 
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Application pt the transformation D results in 


R S, 0c> = R s . , 00 = 


tt O' s in S c - in S c 


where 


S c = S A (t>® S A (t -k) 


for k ^ 0. 

Application of property C results in 


S c =S A (t)eS A (t-k)=S A (t-t 1 ) 


Thus S c is the sequences S A shifted by some delay 1 1 ; by property A,L = 2" -1 
and by property B ' 


ft 0’s - # l*s in S A ( ) = - 1 


and R (k) = - 1/ (2 n -> 1) for k ^ 0; for k = 0 


’ S c = S A (t) ® S A (t) = 2 S A |(t) = 0 mod 2 


thus S c is the all zero sequence (# of 0 f s = L and # l f s = 0) and 


R* (0) = — = 1 ! 
s a v / L 1 


as expected. 
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Cross Correlation of L inear Maximal PN Sequences > : r | 

The cross correlation of two distinct linear maximal PN sequences S A and 
S B is given by equation (1) with 


S c = S A (t) ® S B (t - k) any k 


Unfortunately the shift and add property cannot be used here. Thus we cannot 
predict at this point the values taken by (k) as a function of k. The literature 
is not helpful either in this respect; however some bounds* have been given. 
These bounds are very loose (i.e., the observed autocorrelation values are much 
less than the bound) except for a specific one derived by Gold [23 and applicable 
to only certain pairs of sequences. The results of this reference will now be 
stated and explained. 

* 

Let P (x) and P (x) be two primitive polynomials of degree m with roots 
a and fi respectively and let T (an integer number) be — 



RabOo* 1 /*: 2 "- 1 ) 


This is the only tight cross correlation bound that is available at the present and, 
at the risk of being repetitious, it should be emphasized that this bound applies 

* These are upperbounds on R$ a Sq (It) and imply that this R is less than a certain value (the. 
bound). 
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only to a few chosen pairs of PN sequences. As an example consider m = 5, the 
number of maximal sequences is (from table 1) 6, the length is 31, and T = 9. 

If sequence S B is chosen with generating polynomial (x) of root/3 , then choose 
sequence S A corresponding to polynomial I> (x) whose root is a = ft ; the cor- 
relation (k) for any k will be equal to or less than 9/31 = .29. Notice that 
the roots of (x)* area , a 2 , a 4 , a 8 , a 16 , the roots of (x) are /3 = ct? , /3 2 = 
a i8 t ft = a s , ft a 10 , /3 16 = a 20 (five roots for each polynomial of degree 5). 
There are four more maximal sequences from the four polynomials 


P c (*>* Pj}( x )» P E (*). Pp( x ) 


the roots of P c 

(x) 

are 

7 = 

a 3 , a*. a 12 , a 24 , a 17 

V 

the roots of P D 

(x) 

are 

8 = 

\ 

a 7 , a 14 , a 28 , a 25 , a 19 

the roots of 

(x) 

are 

e = 

a 11 , a 22 , a 13 , a 26 , a 21 

the roots of IJ, 

(x) 

are 

c o = 

alS, a 30 t a 29 t a 27 ( a 23 

Also notice that 


y T : 

= (a 3 ) T = a 27 


thus (x) and (x) are also a pair of polynomials whose sequences exhibit 
the Gold bound. In a similar fashion we can find the remaining pairs which are 
P b *d * p d ^ ^ and^ and finally P p and % . 

Thusforanyk 

*DA( k ) 5 9/31 • • ‘ V ' , 

RpcCk) 5 9/31 

’ : v ’ ; "' < 9/31 • ' 

. _ . ; R^k) 19/31 

V k > £ 9/31 

" ii ■ > i — ' , * ;'i* i'" 

* See Appendix B. . i .. . .. 
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There are two unsolved problems here. The first problem is: What is the 
maximum crosscorrelation of the remaining 9 pairs such as the pairs (A,C), .. 
(A,E), (A,F), <B,C), (B,E), (B, F), (C,D), (D,E), <P,F). The second problem is .. 
what is the distribution of R (k)? 

AB - V — / 

/ 

Presently the answers to the aboye problem, are obtained by brute exhaustive 
simulation. Some shortcuts can be taken by using some of the theory, but the 
shortcuts simply mean less computer time. An exhaustive simulation was per<- • 
formed by the author with m - 7 and 11; the results for m = 7 are outlined in 
the next section,* 

Experimental Results on the Crosscorrelation of Linear Maximal PN Sequences 
of m - 7 


The eighteen PN sequences with m = 7 and L - 127 have been generated f 
and the crosscorrelation of each pair was found exhausively for all k (0 to 126), 
The results showed that for a certain group of pairs the crosscorrelation takes 
on a certain number of values with the same frequency. For example, for a 
certain group * which I shall call Group A the crosscorrelation took on the values 
-17 for 28 k’s, 15 for 36 k's and -1 for 63 k's. The different groups and the 
cross correlation values have been tabulated in Table 3. It can be seen from 
Table 3 that there are three sets of cross correlation values. For certain pairs 
of sequences the crosscorrelation takes on the values -17, -1, 15 with the fre- 
quencies shown under Group A, for other pairs it takes on the values -41, -17, 

-9, -1, 7, 15, 23 with the frequencies shown under Group B and for pairs con- 
sisting of any sequence and its reverse the crosscorrelation taken on the values 
indicated by Group C. 

Table 4 shows how 16 of the 18 sequences divide into Groups A and B when 
each of the eighteen sequences is the sequence being crosscorrelated (its reverse 
belongs to Group C). An explanation of this division into the two groups is given 
in Appendix C. 

§ 

It can be shown that the grouping of the sequences as observed in Table 4 
is strictly a function of the roots of the sequence generator polynomials* . 


1 See Toble 2 for identification of each sequence and its roots 

t This group is the one that includes the a 1 1 polynomial. 

* The results for m =9, 10, 11, and 12 will be published soon. 

S See Appendix C 

♦ The roots of the polynomials generating the sequence are referred to here 
as the roots of the sequence for the sake of brevity, 

, \ 
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Table 2 


The 18 Maximal sequences with m = 7 and the corresponding roots 


Sequence number 
and corresponding 
root of generating 
polynomial 

The complete set of roots 
satisfying the sequence 
(only exponents of a are shown) 

1 - a 

1, 2, '4, 8, 16, 22, 64 - 

2 - a 63 ' / 

63, 126, 125, 123, 119, 111, 95 

3 - a 3 

' 3, 6, 12, 24,48,96, 65 

4 - (a 3 ) 63 = a 31 

31, 62, 124, 121, 115, 103, 79 

5 - a 5 

\ 5, 10, 20, 40, 80, 33, 66 

6 - (a 5 ) 63 = a 61 

\ 47, 94, 61, 122, 117, 107, 87 

7 - a 7 

j ■ 

j 7, 14; 28, 56, 112, 97, 67 . 

8 - (a 7 ) 63 = of 5 

: \ 15, 30, 60, 120, 113, 99, 71 

9 - a 9 

9, 18, : 36, 72, 17, 34, 68 

10 - (a 9 ) 83 = 

55, 110, 92, 59, 118, 109, 91 


11, 22, 44, 88, 49, 98, 69 

— llT- (a 11 ) 63 = a 58 

39, 78, 29, 58, 116, 105, 83 

13 - a 13 

' 13, 26," 52, 104, 81, 35, 70 

14 - (a 13 ) 63 = a 57 

23, 46, 92, 57, 114, 101, 75 

\ 15 - a 19 

19, 38, 76, 25, 50, 100, 73 

16 - (a 19 ) 63 = a 54 

27, 54, 108, 89, 51, 102, 77 

17 - a « 

21, 42, 84, 41, 82, 37, 74 

18 - (a 21 ) 63 = a 53 

43, 86, 45, 90, 53, 106, 85 

Sequences numbered with odd numbers are taken from [3] ; sequences 

whose number is even and next to an odd number are the correspond- 

ing reverse sequences, i.e., 2 is the reverse of 1, 12 is the reverse of 

| 11, etc. 
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/' Table 3 












v k i 

Specifically if the root of a sequencers 77 then the sequences with v belong, to 

group A, the sequences with roots 77 ^ belong to Group B and the sequences with 
roots 77 , ^c belong to Group C. ...... 


For 


. m = 7 (L = 127) 


l = 3, 5, 15, 9, 11, 39, 13, 23, 27, 43 

A 

■ - /. ■ , 

/ 

4 =31, 47, 7, 55, 19, 21 

B 

and ' i 

<e c = 63 x 


Note that for this case 9 n* 17 mod 127 represents the Gold sequence, and 63 
represents the reverse sequence. / 

It should be observed that the crosscorrelation is always an odd number. 
The reason is the fact that L (= 127) is odd, and the # of zeros plus the number 
of 1 f s must equal L while 


LR. = Q . = # of O’s - # l’s = 2(# 0’s) -L odd QED 


Another observation is /chat the cross correlation in group C takes on values 
equal to 3 mod 4. This indicates that the number of zero's in the resulting se- 
quence, S c = S A ® ^ is always odd. (The crosscorrelation in group A is 3 
mod 16 and for grdup B it is 3 mod 8 .) * 

The Correlation as a Random Variable 

/ ■ 

In view of the results of Table 3 the crosscorrelation of the three groups 
will be treated as a random variable. The reason for doing this is the following. 
Although Group A has the smallest maximum crosscorrelation Groups B and C 
seem to be more toward the center. First let us find the mean and variance of 
the autocorrelation if the autocorrelation is regarded as a random variable. 


*See Reference [4] . 
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Table 4 , . 

Correlation Groups for the 18 Sequences of m = 7 ' 

Part 1. If S a is sequence #1 (reverse sequence is' #2) 

Group A = Sequence #’s 3, 5, 8, 9,- 11, 12, 13, 14, 16, 18 
Group B = Sequence #’s 4, 6, 7, 10, 15, 17 

Part 2* If S fl is sequence #2 (reverse sequence is #1) 

Group A = Sequence #'s 4, 6, 7, 10, 11, 12, 13, 14, 15, 17 
Group B = Sequence #'s 3, 5, 8, 9, 16, 18 

Part 3. If S is sequence #3 (reverse sequence is #4) 

Group A = 5, 6, 8, 9, 11, 12, 13, 16, 18, 1 
Group B = 7, 10, 13, 15, 17, 2 , . 

Part 4. If S is sequence #4 (reverse sequence is #3) 

Group A = 5, 6, 7, 10, 11, 12, 14, 15, 17, 2 
Group B - 8, 9, 13, 16, 18, 1 


Part 5. 


If S is sequence #5 (reverse sequence is #6) 
Group A = 8, 9, 10, 11, 14, 15, 18, 1, 3, 4 ■ •" 
Group B = 7, 12, 13, 16, 17, 2 - • ' 



Part 6. 


Part 7. 


Part 8. 


Part 9. 


If S a is sequence #6 (reverse sequence is #5) . 

Group A = 7, 9, 10, 12, 13, 16, 17, 2, 3, 4 

Group B = 8, 11, 14, 15, 18, 1, ' / ; * 

If S is sequence #7 (reverse sequence is #8) 

Group A = 9, 8, 10, 12, 13, 15, 16, 17, 2* 4, 6 
Group B = 11, 14, 18, 1, 3, 5 f 


If S a is sequence #8 (reverse sequence is #7) 
Group A = 9, 10, 11, 14, 15, 16, 18, 1, 3, 5, 7 
Group B = 12, 13, 17, 2, 4, 6 


If S a is sequence #9 (reverse sequence is #10) ; ■ 

Group A = 12, 13, 16, 18, 1, 3, 5, 6, 7, 8, 

Group B = 11, 14, 15, 17, 2, 4 fti l 


Part 10. If S a is sequence #10 (reverse sequence is #9) .. ir ,.j. rn .y- ! 
Group A * 11, 14, 15, 17,2, 4, ‘5, 

-...(Group B= 12, 13, 16, .18, 1, 3 .. .. v;;: , ...... Vo 

Part lli If S is sequence #11 (reverse sequence is ; #12) i ’ iv' 

:• u-c'.'.' i Group A = 14, 15, 18, 1, 2, 3, 4, 5, 8, S. c-wni; i^ocfu/. 

■” Group B = 13, 16, 17, 6, 7, 9 ■' -i »>ro xi a r:;:m 

> • . < ■ iv !>..\ ..i: 

Part 12. If S a is sequence #12 (reverse sequence is #11) 

Group A = 13, 16, 17, 1, 2, 3, 4, 6, 7, 9 - _ 

Group B = 14, 15, 18, 5, 8, 10 !>! 
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Table 4 (continued) 

^Eart"l3. If S is sequence #13 (reverse sequence is #14) 
Group A = 16, 17, 18, 1, 2, 3, 6, 7, 9, 12 
Group B = 15, 4, 5, 8, 10, 11 - 

Part 14, If S a is sequence #14 (reverse sequence is #13) 
Group A = 15, 17, 18, 1, 2, 4, 5, 8, 10, 11 
Group B = 16, 3, 6, 7, 9, 12 

Part 15. If S Q is sequence #15 (reverse sequence is #16) 
Group A. = 17, 18, 2, 4, 5, 7, 8, 10, 11, 14 
Group B. = 1, 3, 6, 9, 12, 13 

Part 16. If S is sequence #16 (reverse sequence is' #15) 
Group A. = 17, 18, 1, 3, 6, 7, 8, 9, 12, 13 
Group B. = 2, 4, 5, 10, 11, 14 

Part 17. If S a is sequence #17 (reverse sequence is #18) 
Group A = 2, 4, 6, 7, 10, 12, 13, 14, 15, 16 
Group B = 1, 3, 5, 8, 9, 11 

Part 18. If S is sequence #18 (reverse sequence is #17) 
Group A = 1, 3, 5, 8, 9, 11, 13, 14, 15, 16 
Group B = 2, 4, 6, 7, 10, 12 



K =: 2/(2" - 1) = l/(2 m - l) 2 

f(0.) is the frequency of the i th value of the random variable. The variance of 
8 becomes: , . , . i . .. . 
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<j\ = 


“v 1 

2_«? 


/(2 m - i) - e 




= (2™ - l) 2 + (- l) 2 (2 m -2) -& 2 _ (2 m - l) 2 + 2 m - 2 - 1/(2"? - 1) 


2 m - 1 


2 m - 1 ' 


2 = <r| /(2“ -l) 2 = 


■ v 1 


2 _ ( 2 m - l) 2 + 2 m - 2 - 1/(2“ - 1) . - ' 


(2® - l) 3 


for the case m = 7, 2“ -1 = 127 and 


V - 1/127 


and 


and 


= j A 12 ?)** 


126 - 1/127 


A , 

: 'H'h -• i 

. ’ vy/f/'r) 


' 

’ y" 

•• ,ru--' 


127 


°r = = -0891 = .00793 

R 127 R 


= 11 . 313358 , 


; . . * i »f v. • ,’C ■; U *' 


Similarly a sjmple calculation reveals that 


\ 




0< A > = 1/127, |0< A > I = 1079/127, = 8.496062, cr£ A > = 11.3133 


< B > = 1/127, i^ B >! = 1133/127 = 8.92, cr£ B >= 11.3133 or cr 6 


_ j/ 16255 


127 


0< c > =1/127, |0< c > | = 1217/127 = 9.5826, <r$ c > = 11.3133 

•' > ■ ’ . ■’ ’ y : 'V. :.i. •< 

where the superscripts A, B, C t refer to the groups A, B, C respectively. 
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The interesting result is that the mean and variance of the autocorrelation or i 
crosscorrelation for groups A, B, and C are independent of the group. The mean 
of the absolute value is approximately 



less than 10% 


This value is not bad when' one is not concerned about peak crosscorrelation, 
Since R*is the same foj / Groups A, B, or C and since <r R is the same also fc 
Groups A, B, C, comparison of A, B, and C can only be done on basis of _|R| 
peak and [r|. On both 7 these criteria group A becomes superior to the other 
groups. We call G/oup A a group with good crosscorrelation properties. 


Sets of Sequences with Good Crosscorrelation Properties. 


From Table 4 one could choose a set of 3 sequences with good cross - 
correlation properties, i.e., all R R R belong to a set where the pair- 


wise crosscorrelations are three valued as shown in Table 3 under Group A. 


For four signals one could similarly choose 1, 3, 5, 8 with R 13 R 1S R 18 R 3S 
R 38 and R sg all belonging to a set having the same crosscorrelation distribution. 

For 5 signals one could choose 1, 3, 5, 8, 9. For 6 signals 1, 3, 5, 8, 9, 11 
breaks down because (9, 11) belongs to Group B, but 1, 3, 5, 8, 9, 18 works. 

For 7 signals one must consider many combinations but in general one can 
be sure that no set can be found of more than 9 signals such that for all i, j. 
i, j = 1, 2, . . . , ten or more 0. . belong to Group A. 

On the generation of a large set (2" 1 + 1) of non-maxi mal PN sequences with 
good 'crosscorrelation properties. The crosscorrelation properties of the 
previously described sequences are not at all discouraging even for Groups B 
and C because the only difference is in the peak values which come with a very 
small probability, i.e., for Group B Q - -41 occurs only once out of 127 possible 
times. Needless to say that the autocorrelation of each of the sequences is al- 
most perfect (orthogonal), as shown earlier. 


The only difficulty arises when one needs more than <£(L)/m sequences. One 
solution to this problem is to generate sequences S. from two maximal sequences 
S A and S B of the same m as illustrated in Figure 2. 

V . • .1 ... .. . . . 


* ® = 127 ’ °R = /12 ^* 
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This configuration will generate* 2 m + 1 non-maximal sequences of length 
L = 2 m - 1 whose autocorrelation and crosscorrelation for any two sequences will 
have the identical values (with identical frequencies) as the crosscorrelation of 
S A and S B . In other words by choosing S A and S 0 so that belongs to Group 
A, then S. and S . as generated by Figure 2 will be such that R . . will also be- 
long to Group A lor all i - 1, 2, ..... 2'" + 1 and j = 1, 2, . . . „ J . 2 m + 1* 


The Correlation of PN Sequences Without limitation on the Length. 

It was proven in Part 3 that if the correlation is a random variable for the 
case of a linear maximal PN sequence, the autocorrelation, had a mean equal 
to (1/L) 2 and a variance - 

X. •• , 

I(2 m - l) 2 + 2 m - 2 - 1/(2” - 1) ... • 

. .... • ... (2™ - l) 3 . ’-i-W.-VA . 


which for large m can be ^proximated to 



where L is tht/length of the sequence. 


1 ■ - n V/ 1 ; , 


•:: r i l 

>r Hr i*. i'jrif 

S: ? 0* ' 1 \ • , h J : 


* See Appendix D * : 

** One exception is the autocorrelation of 2 of these sequences $| = ©0 and 

$j =S fl +0 =S B which is — 1/(2 m - 1) for r?0 and 1 for r=0. 
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It has also been observed that the crosscorrelation of any two maximal PN 
sequences of the same degree (or length) has the same variance and mean as 
above. 

The question that arises is: "What is the mean and variance of the auto- 
correlation or crosscorrelation when the sequence is not cr03scorrelated for its 
entire length?" 

a. Letting the segments be p bits along with p < m, the correlation of two p 
bit segments of a sequence S is 


* <t) m£ v°> v t ) 

\i B ° -\ 

\ v 

\ 

\ 

by applying properties D and C 


9 (r) = [# of 0’s - # of l’s] in S Q 

J ’ 

where 

S C = S (0)®S(T) = S (Tj) - 
letting w denote the weight* of the p-tuplet then 


6 (t) = p - 2 w 


Since w is not known for each r, it is treated as a random variable; once the 
probability density** function of w is found one can find 6 and 6 2 . For each w, 
there are C P V p-tuplets since there are ways of having w "ones" in p distinct 
(ordered) positions. 

-- ' .ii ,.i . 

•The weight is equol to the number of l’s - 

** Since w is discrete, insteod of the probability density we find the frequency of occurrence f 

w 
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Obviously the frequency of a given distinct p-tuplet* in a maximal sequence 



Then the expression for 6 becomes: 




P 

mp -cE 


w 2 n-p C p 


- P - 

From Appendix E 


2 m + 1 ~p 


L 


w c; 


£ 


w C p = E 2 P 
v * 2 


and 


consequently . . 


2 HI + 1 *“ p 

T 2 


6 = p - £ 2 p =-E [2° - 1 - 2“] = -- H 


R S £ s .i 

P L 


*For a Shift Register generated sequence the ail zero p-tuplet has a frequency of 2 m " p -l be* 
cause the oil zerom-tuplet is the only m-tuplet which is absent from the maxima! sequence. 
However, in this discussion this fact will be neglected for the sake of simplicity with no 
considerable error. 

** This is the amount of times that a distinct p-tuplet appears in an nvdimensionai space over 

GF > 
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The standard deviation of 6 t <?% , is be definitions 



the first term becomes 


now 


2 m ~ p 
L 

p 

L 








i w = (i + i) p 


= 2 P 




from Appendix E 




The first term becomes 


2 ®“? 





2 P - 2 • 2 P +- 
P 


2 ? 



= 2 m /L 
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If p is greater than m (with still linear sequences) these relations, although 
no longer true, may be considered as upper bounds for the following reasons. 

If p - m = b, given the first m bits the remaining b bits of the p-tuplet are 
determined depending on the feedback connections of the shift register; further- 
more there are only 2 m p-tuplets (not 2 P ) which obviously makes some* p- 
tuplets inadmissible. The majority of the above inadmissible p-tuplets are the 
type whose weights is much different from the mean weight w ; (i.e. p-tuplets con- 
taining more than m consecutive 0*s or l's). Consequently the standard deviation 
should be smaller than the one found above. 

Now if one assumes that the above inadmissible p-tuplets exist in the sequence 
but that their frequency is 2 m ~ p (which in this case is a fraction of unity) the 
equations of section 5 should still hold so far as the frequency of a distinct 
ordered p-tuplet is concerned, (These equations should not be considered an 
upper-bound) i.e. if only one p-tuplet out of the 2s*" n (p > n) is admissible we 
may say that on the average the frequency of a distinct p-tuplet is 2*“p or 
1 / 2 p -m (p > n ) go th a t the frequency of a given w is still 2 n ~ p C p . . 


* 2 p * , m 
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APPENDIX A 


This appendix is used to arrive at the formulation of the transition matrix . 
A by treating the contents of the shift register as a state vector. The analysis 
has many gaps. For more rigorous analysis the reader is referred to [1]. 


Definition: A maximal pseudo-random sequence S A is a sequence generated 
by a feedback shift register connected to divide by the primitive polynomial P A 
(x); the period is 2 n - 1. P A (x) is of the form . * 


P A (x)= £ h,** 

i - 0 


with h. = 0, i and h Q = h n =1 


( 1 ) 


and it is the characteristic equation of -the transition matrix A that gives the re- 
lationship between the state (contents) of the Shift register generator at some 
time t and another time t +r . If u(t) represents the contents of the shift register 
at time t then . . ' 


( 2 ) 


(3) 


u(t) is an n-dimensional vector and if it is considered as an n-symbol segment 
of the sequence S A then u(t + 1) will contain n - 1 of the symbols in u(t) and a 
new symbol which is a function of the n - 1 previous symbols as shown by the 
matrix A in the companion form. ; 


u(t + r) = A T u(t) 




or 


u(t + 1) = A u(t) 




0100000 . . 
00100 . . . . 
0 0 0 1 0 . . . . 


0 0 
0 0 
0 0 


0 0 0 0 .. .... 0 0 0 1 I I u n- 1 (t)| 

h^h^. . . . VA-»I |u n (t) 

The above equation results in the following relationships 


V 1 ) 


, ( 3 a) 


Wi (t + 1 ) = u. +l (t) t — 1 , 2 , . • - • n 1 


and 


\ 


n - | n » I 

u „( t + 4) = 2Z hiUi+iCt) = YL h i u i (t + ^ 


i = 0 


i =0 


( 4 ) 


( 5 ) 


Since the last equation 


n - 1 


u n ( t + 1) = 2^ w u i ( t * j > 

i =r O 

:\ 


( 6 ) 


is independent of time we can drop the time dependent thus resulting in the re- 
cursive formula 
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( 7 ) 


i 



U n = 



i “ 0 


h. U r 






which gives the nth symbol of the sequence in terms of the n previous symbols. 
It should be noticed that the h L *s of equation (1) are identical to those of equa- 
tion (7). ; j 


\ 
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APPENDIX B 


ON THE ROOTS OF POLYNOMIALS OF A SPECIFIED PERIOD 

This appendix is written in order to prove some interrelationships between 
roots of polynomials. 

1. Let ' ‘ - 





+ 0 = P(a 2 > 


because 


h 2 = h A in GF (2) and 2 = 0 mod 2 
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I 


Thus 


t ; < 


: . . t 4 


P (a 2 ) = fP(a)] 2 = 0 




We conclude that /3 = a 2 is also a root of P (x). In general if k = 2 


[P(a)] k =P(a k ) = 0 


because 


but 


hj a* = h j ail{ + C j (other terms)-. 


I— 1 


, k ! 

C, = T~ T-" t = 0 mod 2 for j * k. 

1 (k-j)!j! 


V' 


In general if a polynomial has a root a then a 2 is also a root» When T' > 2 m - 1 
then 2^ is taken mod 2 m - 1 which is the period, of the root. Example for m = 5 
2* - 1 = 31. If jB = # W is a root then are also roots of- the polynomial; then 
the roots are: > 


a 


li 


22 a 44 mod 3 6 


~ a. 13 . 


,26 


and 


a S2 mod 31 = a 21 


thus any of a 11 , a 22 , a 13 t a 26 or a 21 satisfy the same polynomial of degree 5. 

2. The Roots of the Reverse Polynomial. 

The primitive characteristic polynomials can be found in Tables. Some of 
the tables contain only half of the primitive polynomials. The other half are the 
ones representing the reverse sequences and can be obtained by finding the 
polynomials due to /3 201-1 / l for each fi . Proof: Let ft be the root of an arbi- 
trary primitive polynomial* 

/ / i 

/ "■ P(x) = X m + X p + X q + X +1 
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then 


P (ft) = ft u + p* + 


/3 q 


+ l_a 0 



Now consider the polynomial due to the reverse sequence* 


P R (x) = x u + x n_p + x n “ q 




1 

We shall test the root \ 

' \ 



ft2 a ~ 1 “ 1 - 0-2*' 'mod 2° ” l 


on P R <x) and more generally we shall test the root ft' 2 on P R (x) k = 0, 1, 2, 

m - 1 because / _1 , ft' 2 * ft'l ft' 2 \ . . . . ft * 3 "* 1 are roots of the 

same polynomial {ft' 1 - ft “ 2 '’’ n,od 2 ° ~ 1 = /3" 1 ). Thus the polynomial P R is 
satisfied by ft' 1 ; it is also satisfied by the above mentioned roots. 


J 1 1 ' 

P R (/T*> = - + — + — + 1 
v, v ft n ft*~p ft n ~<i 

.* , . 

= ft' n [1 + /3 P + ftP + ftt + ft n ] 


but the equation in the brackets is zero (P(/3)). Thus 


P R ( ft" 1 ) = ft~ a P</5) = 0 QED. 



*Sea [ 1 ] for proof. 
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t 


k 1 

Consequently /3 2 is a root of P R (x) \ 

r ■. 

for k = 0, 1, m - X. Thus the positive exponents of the roots that satisfy 

P„(x) are t . . \ - ■ , ■ . \ < • 


2 m -;1 - 2 k , . k.= .0». 1,2 -r. -vii 
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APPENDIX C } 

CROSSCORRELATION GROUPS AND ROOTS OF SEQUENCES 

This appendix explains the grouping of sequences^ for S a equal to one of 
the 18 sequences. The first part of Table 4 gives the sequences in Group A and 
in Group B when S a is sequence #1 (its reverse is #2) and it is rewritten for 
convenience. 

S a = Sequence #1 (root is a ) 

GROUP A = Sequence #’s 3, 5, 8, 9, 11, 12, 13, 14, 16, 18. 

GROUP B = Sequence # f s 4, 6, 7, 10, 15, 17 

Another way of identifying the sequences is by using the exponent of one of their 
roots. Since each sequence has 7 roots we shall use the smallest of the 7 ex- 
ponents (the first entry of the second column of Table 2). Thus the sequence 
number and root identification become the pairs (sequence #, root exponent) = 
(1,1), (2,63), (3,3), (4,31), (5,5) (6, 47), (7,7), (8,15), (9,9), (10,55), (11,11), (12,39), 
(13,13), (14,23), (15,19), (16,27), (17,21), (18,43). 

On the root identification basis, when the root of S a is a (root exponent = 1) 
Group A becomes the sequences whose roots have exponents 3, 5, 15, 9, 11, 39, 
13, 23, 27, 43 and similarly Group B becomes \ 


B = 31, 47, 7, 59, 19, |21 


Group C consists of the reverse sequence only with root* a 63 . Notice that the 
sequence whose roots is a 17 = a 9 mod 127 is in group A as predicted by Gold t2i 
Notice also that this is not the only sequence whose crosscorrelation with se- 
quence #1 is less than or equal to 17. Nine more sequences have identical cross- 
correlation values and frequencies. I 

f 

Now consider the second entry, i.e. sequence #2 with root identification 63. 
Setting a 63 = 0 the root of sequence #2 let us find/3 3 , 0 s , 0 1S , 0 9 , 0 11 , 0 39 , 
/3 13 ,/3 23 , ^2 7 and /3 43 . Notice that the exponents of the above roots are the ex- 
ponents of the a’s in.case 1 for Group A. 
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fi 3 -(a 63 ) 3 '- a 31 


sequence # 4 


< <. .. * • f 


P 5 = (a 63 ) 5 = a 61 < *■ sequence # 6 


/?»S 

/6 9 

/8“ 

/J13 

/323 

/3 27 


= a 56 < > sequence # .7 

_ a s 9 ^ > sequence # 10 has root y3 17 

= a 58 < ► sequence : # 12 


= a 11 *e- 


sequence #11 


= a 57 <— — > sequence # 14 


- rt 52 


•= > sequence #15 J 

// . - " 

= a. 42 < * sequence # 17 


*■ sequence # 17 


but sequences #4/6, 7, 10, 12, 11, 14, 17, 15, 18 form Group A for sequence 2 as 
seen in Table 2/ This result shows a possibility that the crosscorrelation Group 
of values have/a dependence on the roots of the sequences. Thus if a is a root of 
the sequence/S a then Group A consists of sequences whose roots are a^ A with 
A given bj 

X ■ ' • • • 

= 3, 5, 15. 9. 11, 39, 13, 23, 27, 43 

k ’ ■ / - • ‘ * • 4 • - ■ • 'ir* 7 /=•:•■* ’ ‘ 

- . ■ 

and group B consists of sequences whose roots are o' 3 with l given by i: - •' 


-e B = 31, 47, 7, 55, 19, 21, ; V/..4/ 


t>V ( * 


i'f .’ ; 1 : ‘ ! x 


similarly ’ ’ ’ P >• *•-** 

•/ ■ ; ii.- K;\> :■.!*' r-~ {•■ , yi.' 1 '.’ •/.;•!; y; 1 ‘ ' 

. -Vi* 1* 

'A = 63 
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As we have seen if is sequence #2 with root /3 ( /3 = a 63 ) then to find group A 
all we need do is raise /3 to the powers l A as we have done, with results identical 
to the ones in Table 2, part 2, Group A. Thus for Group B raise /? to i h and 
identify the corresponding sequences, , 


/3 31 

= (CL 63 ) 31 

n 

a 48 



sequence 

<# 

3 

/3 47 

= (a 63 ) 47 

j 

a 40 

«- 

- 

sequence 

# 

S 

fi 7 

= (a 63 ) 7 


a ls 

— 

- 

sequence 

# 

8 

/3 SS 

_ (a 63 ) ss 

= 

a 36 

— 


sequence 

# 

9 

p i« 

= (a 63 ) 19 

= 

a 34 

- 


sequence 

# 

16 

£21 

= (a 63 ) 21 

zz 

a 53 



sequence^ 

# 

18 


K' " . 

• f \~ 

The sequences are #'s 3, 5, 8, 9, 16, and 18 and they agree with the Group B 
entry in part 2 of Table 2. (Of course Group C is sequence 1; check: 


£63 = (a 63 ) 63 = a 32 * *■ Sequence ft 1). 


Let us derive now the Correlation Groups for the 6th sequence. Let the 
root of sequence #6 be 7 - a. 47 Table C gives the information on the powers of 
y , a , and corresponding sequence number for Groups A, B, and C. It can be 
seen that the numbers at the third column correspond to the ones for Groups A 
and B in Table 4 and Group C is sequence 5 as expected. The remaining entries 
of Table 4 could be derived in this manner. v 
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TABLE C 


y - a 47 <s >• sequence #6 



Exponent of y 

Exponent of a 

' Sequence # 


3 

3(47) = 14 mod 127 

7 


5 

5(47) = 108 mod 127 

16 


• 15 

/•' 

15(47) * 70 mod 127 

13 


/ 

9 , . 

9(47) = 42 mod 127 

17 


11 

= 9 mod 127 

9 

Group A, 

39 

= 55 mod 127 . 

10 


: ” 13 ".' 

= 103 mod 127 

4 


23 .... 

. = 65 mod 127 

3 

X 


27 

/ = 126 mod 127 

\ 

= 116 mod 127 

■x/ 2 

> r \ } 4 * ,■ 

43 

r j j . t v . . . . . ; t , , 

12 

i # 

31 

= 60 mod 127 

8 


47 

= 50 mod 127 - 

15 


; ‘ 7 : ' 

"» 75 mod 127" 

14 

Group B 

59 

'v 

. , • 19 .,,•/• ...... 

= 106 mod 127 
■= 4 mod 127 

’ 18 

. ... ... 1 • :. 

/ • '* ; ; : 

Z1< ■ ■ 

= 98 mod 127 

; 11 1 • " 

Group C , 

j,’ ’ 63 ‘ .... " 

= 40 mod 127 

5 ... • 


• >i ; : ‘ t . 1 {.■ - : 1 • 
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APPENDIX D 


This appendix is devoted to the development of certain sequences S. of Fig. 

D1 with identical crosscorrelation properties as the original two generator se- 
quences, Before proceeding we repeat four important properties that will be 
used in the course of the proofs. j 

| 

A Sequence of period L can be viewed as a sequence of L states (n-dimensional 
vectors L = 2 n - 1 for binary). j 

r 

Property 1. A maximal binary sequence will occupy all 2 n - 1 non-zero states 
before it starts repeating. j 

i 

S 

Property 2. (Balance property). The number of l’s in a maximal sequence 
exceeds the number of 0’s by 1. j 


Thus the number of zeros of a maximal sequence is equal to 2 n ~ 1 - 1 and the 
number of l's is equal to-2 n ~ K . • • - 


Property 3. (The Shift -and-add property). 




When a maximal sequence is first shifted by t bits and then added to itself 
the resulting sequence is the same sequence delayed by some time t i.e. 


S A (°) + s A (t) = s A ( tl ) t /0 


* 


This property becomes obvious from the following argument. A shift register 
generated sequence described by equation (7) in Appendix A in a linear 
process. Thus if the elements of S A (0), u n obey the equation 


u n (0> = 7} h i u i<°> 

i *0 

then the elements of S A (t) obey the equation 

\ 


*AII additions ore mod 2 unless they result in some statistic R or $ ond their functions. 


O 
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/ 


n - 1 

«»(*> = h i u i(t) 

/ 1 


i = 0 


thus the elements of the sequence [S, (0) + S. (t)] are 

ri A 




n - I 


V n ? u n(°) + u n ( t ) = h i ^ U i C°) + u i( t >^ = h i V ' 


where \ 


i =0 


i a 0 


V. = Uj (0) + u 4 (t) 


Since v n obeys the same equation as u n (the same set of h.'s) the sequence whose 
element is v n is the same as S A except for a delay, if S A is maximal. (If S A is 
non-maximal then v n may be another non-maximal sequence obeying - ;he same 
recursive equations). 


A transformation that changes the correlation function to the difference of 
the numbers of l*s and zeros in a sum sequence. Consider two sequences S A and 
S B . The correlation of S A and S B as a function of t is defined as 


defining 


v 

' ' ' i I i 


S,(0)-S E (t) 


«»,(<)= LI B (t) 


AB V 


(t) turns out to be the inner product of S A and ^ (t). *11103 


Thus 


**,<0 = 0 tS A (°>, 


/ V’ ■ ■ i ' ; i I j ( .**[/; 


* 2_ “■ • b H. 


V . ’ ' it. 
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where 


— [&ji flji 0 ^ • « • * 1 1 Sjj — (b|t bj* b 3 


. . ] and aj, bj = - 1» 1 


if a t , b j =1,0 then one can apply the following transformation 

1 -* 0 - 1 -• 1 resultingin * -* ® ( mod 2addi t ion) 


(Another way to look at this is in terms of agreements and disagreements). The 
summation sign adds all agreements (l’s) and subtracts the disagreements (-l’s); 
now that we transformed the l’s to 0’s and the -l’s to l's, we must add the zeros 
and subtract the ones of the sequence resulting from the modulo 2 addition of 
S A (0) and Sg (t). The reason for this transformation is the fact that very often 
we are able to determine the counts of 0’s and l’s of a sequence resulting from 
a mod 2 addition of two other^sequences, without going through the mechanics of 
adding bit by bit. Subsequently this count of 0’s and l's will readily enable us to 
determine the correlation of S A and S B (t). If .. 

/ = + S B (t) 

then / * 

^AgC t) - [the number of zeros inS c ] - [the number of ones in 


By using the above properties we shall now determine the autocorrelation function 
<9 m (t) of maximal shift register generated sequences. Letting S A be this sequence 
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where 


# M (t) = (ft of zeros in S c ) - (ft of ones in S c ) 



S c = S A (0) e S A (t) 


by using the shift-and-add property for the maximal sequence S A . , 


S c = S A (t f ) for tV 0 


Furthermore using the balance property for maximal sequences 


resulting in 


• . r , .y 

[# of zeros in S A (t')] - f# of 1's in S A (t'). ] = ■- 1 

• : • ■ :< 

^(t) = - 1 for ty 0 : •' 


or 


X 


R AA (t) = 


- 1 


2 n - 1 


for t = 0 


S c = S A (0) + S A (0) = 2 S A (0) =0 


resulting in the all zero sequence with 2" > 1 zeros; then 


* M (0) = [# of zero in S^J - t# of ones in S c ] = 2 n - 1 - 0 
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and 


Raa(O) = <2" - l>/(2" - 1) = 1 


Suppose now that we generate sequences S by the polynomial P(X) which is the 
product of P A (x) and ^ (x), the primitive polynomials of degree n which were 
defined previously. 


P(x) =P A (x)-P B (x) 


This is equivalent to taking S A (t) and Sg (0) and mod 2 add them (bit by bit) 


Si CO) = S A (tj) • SgCO) 

This can be accomplished by the configuration of Figure 1. 



Figure D 1 


It will now be shown that the auto and crosscorrelation of the sequ ences S { 
generated by the technique of Figure 1 takes on the same values as the cross- 
correlation of S A and S fl . The cross cor relation of two sequences S. and S. 

(i ^ j) is 1 

’ ' \ • ' 

0 tj (t) = tfCS^O) e S^t)) = # of 0’s in ^ - ft of l’s in ^ 
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where 


Sj< = ^( 0 ) ® Sj(t). , 

/ 

Substituting for and S. 


S i ( 0 ) = S A (t i ) * 83 ( 0 ). 



S i<t) =s A (tj + t) ®S B (t) 


l 

i 


- / s k(°) = s A (t.) * s A ct j + t) *^(0) ® S 8 (t) 


= S A (^) e S B (^ 2 ). 


Now going backwards the numbers of zeros and ones in S k is equal to the 
numbersjot-zeros and ones in 

W +S B ^ thUS 


[# O’s in S k ( 0 )] 


.£# of 1 ’s in S k ( 0 >] = Sg C^ 2 )J 


or 




thus the crosscorrelation of S t and S. is equal to the crosscorrelation of S A 
and S B at some other time delays. Let us now examine the upper bound of the 
autocorrelation of 


6 ii (t) = tf tSjCO). S i ( t)] = [#of O’s in S tl ) - 


[ # of l’s in S^] 


and t 0 
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where S . . = (0) + S . (tt,/ By making the proper substitution 


Ki r S A (t.) •+ S B (0) « S A (t. + t) •-S b (t> 


= S A (t A > + S A (t i + t) © S B (0) + S B (t) . t f 0 

■ '• ’v 

= S,(V+ S,(-C 2 ) 

- / > . 

and .. 




Yielding the same results as before. Thus both autocorrelation and cross - 
correlation of the sequences S t generated by the products of the primitive poly 
nomials P A (X) and P B (X) are equal to some crosscorrelations of S A and Sg . 


For the case where S. , the interfering sequence, is really a linear sum of 
other sequences, (but notS. the primary sequence), the crosscorrelation then 
still has the same values as the crosscorrelation of S A and S g . This will be 
shown in the sequel 


let 


x. 




s «e < mod 2 


The sum sequence S v becomes 

& 


but 


Sj, = S. © 


J2 S^(t^) mod 2 

I 



S i = 


S B (°) ® s A 
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and each S^(t^) is the sum of some S b (t^ ) and ) thus 


S ‘ = B A (.) . S B (0) t £ js B . s/ ^j - m6d 2 

' ^- S ° ('•«,) *71 s * (%) ' "* od 2 


Now by using the shift and add property 


V 
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APPENDIX E 


w w 

> t o 

' V tz 0 W = P /2 + 


replace w in the second term by p - k; then the second term becomes 




k-p/2- 


smce 


• k = 0 


O' ■ (, M ■ tnSW ; 




Consequently / becomes 


• \ 


p£2 p/2 

” w (P “ (k) re P laci ng k by w 

kSQ 


w X o 
P/2 


P/2 

<)> 

w = 0 w = 0 


= H'O'L. <p -" ) (») 


P/2 

= (w +CP q = p 2] ( p 

®=0 VBO 


(A-l) 


Now 


£ o* £© a a- .. 

«r a 0 w ■ 0 


(A-2) 
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APPENDIX F 


To find 


p 

L 


i 2 ; letting i 2 = i(i - 1) + i 


the above expression is equal to 


L (D * X l 0 


since the first two terms are 0, (for i = 0, 1) 



let i - 2 ~ k 


p 

L 


• - 2 p from Appendix E 
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p 



2 P + p(P - 1) 

/ 


P - 2 

E 


- (P - 2) ?, 
k ! (p - k - 2) ! 


/ 4 1 ' &-*««*>* 

= \i v +p(p- i) 2 p - 3 * 

/ 

/ 
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